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Abstract 

Radio frequency driven plasma jets are frequently employed as efficient plasma sources for sur- 
face modification and other processes at atmospheric pressure. The radio-frequency driven micro 
atmospheric pressure plasma jet (/itAPPJ) is a particular variant of that concept whose geometry 
allows direct optical access. In this work, the characteristics of the /iAPPJ operated with a helium- 
oxygen mixture and its interaction with a helium environment are studied by numerical simulation. 
The density and temperature of the electrons, as well as the concentration of all reactive species 
are studied both in the jet itself and in its effluent. It is found that the effluent is essentially free 
of charge carriers but contains a substantial amount of activated oxygen (O, O3 and 02(^A)). 
The simulation results are verified by comparison with experimental data. 



I. INTRODUCTION 



Microplasmas operated at atmospheric pressure have recently gained very high attention. 
A series of review papers and topical issues of leading scientific journals have discussed their 
technological promises as well as their scientific challenges [IHH] . A particularly popular class 
of microplasmas are radio-frequency driven plasma jets, originally proposed by Selwyn and 
co-workers in 1998 [6] . Many groups have investigated their prospect for surface modification; 
the studied processes include etching of tungsten, deposition and etching of silicon oxide, 
and cleaning of thermolabile surfaces from contaminants [71-[T^. In order to improve the 
access for optical diagnostics, Schulz-von der Gathen and co-workers modified the originally 
axis-symmetric design of the jet to a plane-parallel geometry |13]. Their device, termed the 
micro atmospheric pressure plasma jet (/xAPPJ), is depicted in Fig. [l] It is a geometrically 
and electrically symmetric capacitively coupled discharge driven via an impedance matching 
network by RF power of nominally 10 W at the frequency of 13.56 MHz. The plane-parallel 
stainless steel electrodes are 30 mm long, 1 mm wide, and separated by a gap of 1 - 2 mm. 
To allow direct optical access, the sides are made of quartz glass. The jet is operated with 
helium as the carrier gas and an admixture of oxygen at a relatively high flow rate of about 
1 slm. The electron density is a few lO^'' cm~'^, the ion densities are up to 10^^ cm^^, and 
the radical densities can reach several 10^^ cm~^. 

This work expands earlier one-dimensional simulation approaches by Waskoenig et al. 
[11]. We study the yuAPPJ by means of a two-dimensional numerical simulation. This 
allows us to study not only the discharge region but also the effluent of the jet which 
expands freely into the ambient atmosphere. In the effluent the gas is characterized by the 
absence of charged particles and a concentration of reactive radicals. It was found that the 
concentration of atomic oxygen decays only slowly outside the plasma volume and that the 
concentration of ozone even increases for a few centimeter. These phenomena have been 
found before by experiments [T5HTS] . (To show the reliability of the simulation results we 
present a comparison with experimental findings of Ellerweg et al. |18j.) 
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II. DESCRIPTION OF THE NUMERICAL CODE 



The simulations of this study are carried out with the fluid dynamics code nonPDPSIM 
designed and realized by Kushner and co-workers. For a detailed description of the code and 
its application to various types of plasmas and gas discharges, see the publications [T9H22] 
and the citations therein. Here, we just briefly discuss the implemented equations and the 
underlying physics. 

The code nonPDPSIM simulates the dynamics of weakly ionized plasmas in the regime 
of medium to high pressure. It takes into account the physics and chemistry of charged 
particles - electrons with mass m^. and charge — e, ions with mass mj and charge qj - and of 
the excited as well as the ground state neutrals (mass mj). For all species j, the continuity 
equations (particle balances) are simultaneously solved, where Tj is the particle flux density 
and Sj the source and loss term, respectively: 

^ = -V-f, + ^,. (1) 

The fluxes are calculated from the momentum balances in the drift-diffusion approximation 
evaluated in the local center-of-mass system. Dj and /xj are the diffusion constant and the 
mobility (if applicable) of the species j, E \s the electrical field, and v is the mass- averaged 
advective velocity of the medium: 

= rijV - DjVnj + ^^fijUjE. (2) 

For the electron fluid, additionally an energy balance equation is solved which takes into 
account Ohmic heating and the energy losses due to elastic and inelastic interaction with 
the neutrals and ions as well as heat conduction, 



dt U 



Hi. (3) 



To capture the non-Maxwellian behavior of the electrons, all electronic transport coefficients 
(the mobility /ic, the diffusion constant and the thermal conductivity Kg) as well as the 
electronic rate coefficients in eqs. ([T| and ([2]) are calculated by the local mean energy method: 
A zero-dimensional Boltzmann equation for the electron energy distribution /(e) and the 
transport and rate coefficients is solved for the locally applicable gas composition and various 
values of the electrical field. The tabulated data are then consulted in dependence of the 
fluid dynamically calculated electron temperature T^. 



The plasma equations are coupled to a modified version of the compressible Navier-Stokes 
equations which solve for the gas density p, the mean velocity v, and the gas temperature 
T. Here, p is the scalar pressure, given by the ideal gas law, and r the viscous stress tensor, 
related to the velocity shear via the viscosity. The contributions to the energy equation from 
Joule heating include only ion contributions; the heat transfer from the electrons is included 
as a coUisional change in the enthalpy. The heat capacity and the thermal conductivity of 
the gas are Cp and k, respectively, Ahj denotes the enthalpy change due to reaction j: 

I = -V-M (4) 

9{pv) ^ _ Y . - _^ ^{qiUi - miHiSi)E, (5) 



dt 

d{pc,T) 
dt 



-V • {-kVT + pvCpT) - p{V-v){t-Vv) - ^ AhjSj + ^ ■ (6) 



Finally, the potential $ is calculated from Poisson's equation. (The code works in the 
electrostatic approximation so that E = — V$.) The charge density on its right stems from 
the charged particles in the plasma domain and from the bound charges ps at the surfaces. 
The coefficient e — So^r represents the permittivity of the medium: 

-^■{eW^)^Ylm + Ps- (7) 

j 

The surface charges are governed by a separate balance equation, where a is the conductivity 
of the solid materials and the subscript s indicates evaluation on the surface: 



dps 
dt 



(8) 



j 

The dynamical equations are complemented by an appropriate set of boundary conditions. 
Electrically, the walls are either powered or grounded. With respect to the particle flow, 
they are either solid, or represent inlets or outlets: The flow is specifled to a given flux, while 
the outlet flow is adjusted to maintain the pressure. Finally, it is worth mentioning that the 
actual implementation of the equations poses some difficulties due to the vast differences 
in the time scales of the dynamics of the plasma and the neutrals. These difficulties are 
overcome by the methods of time-slicing and subcycling. 
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III. THE DETAILS OF THE SIMULATION CASE 



The described model is employed to simulate a /xAPPJ, similar to the one depicted in Fig. 
1. In our case, the electrodes are 30 mm long, their separation is 1.8 mm. The simulation re- 
solves two Cartesian dimensions x and y. In direction z, translational invariance is assumed. 
(This assumption reflects a limitation of the code. We believe that the corresponding mod- 
eling error is not too severe; it is still possible to resolve the plasma distribution between 
the electrodes and the evolution along the direction of the gas flow. The main error of our 
"inflnitely wide" model lies in the boundary conditions of the advective velocity of the fluid 
at the quartz boundaries. It translates into the coefficient of the Hagen-Poiseuille equation 
being off by about 60%, not too severe a problem as we control the total flux through the 
jet, not the pressure differential between and outlet.) 

The simulation domain is depicted and explained in Fig. |2} The triangular unstructured 
mesh consists of approximately 10.000 nodes. 80% are located in the plasma domain includ- 
ing the effluent, their density is highest between the symmetrically powered RF electrodes. 
The outer boundary of the simulation domain is electrically grounded. 

The device is operated as a helium-oxygen /lAPPJ in a pure helium atmosphere at 10^ 
Pa. The inlet is fed with a gas mixture with a high flow rate of 1 slm, corresponding to a 
maximum advective velocity of about 1500 cm/s. The mixing ratio of helium to oxygen is 
1000 to 5; experiments indicate that these conditions maximize the concentration of atomic 
oxygen [13]. To maintain the outer environment, pure helium is injected at a flow rate of 1 
slm at the left side of the simulation domain. The outlet is controlled to maintain a constant 
pressure. The following species are taken into account: Ground state neutrals O2, O, O3, 
and He, 02(2^), representing the flrst four vibrational levels of O2, the electronically excited 
states 02(^A), 02(^S), 0{^D), 0(^S), and B.e{^S), positive ions O^, 0+, and He+, negative 
ions O2 , 0~, and O3 , and electrons. The reaction chemistry applied here is described in 
more detail in [2T]. Following we neglect ionized helium dimers He^ and excimers He2. 
It has been shown that these species (which require a high electron temperature) play a 
minor role when the discharge is dominated by the reactive oxygen species; this holds for 
mixtures with an O2 content of 0.5 % or more. Finally, the secondary electron emission 
coefficient of the surfaces is set equal to 0.1. 
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IV. SIMULATION RESULTS 



As a result of our simulations, we obtain a complete picture of the dynamics of the /lAPPJ 
under certain conditions. In this manuscript we focus on phase-averaged quantities. The 
results correspond to a total RF power deposition in the plasma of 0.5 W. (This value is 
considerably lower than the nominal RF generator power of about 10 W of the corresponding 
experiments. The deviation can be explained by parasitic losses in the matching network 
and the feed cables. Order-of-magnitude differences between the generator power and the 
plasma heating power are also observed in other regimes f^M 1^.) 

Figs. |3] and |4] show the two-dimensional spatial distribution of the electron density ng 
and the electron temperature Tg. The ID profiles (normal to lateral gas flow) of the charged 
particle densities and the electron temperature at a position x=2 cm, are shown in Figs. 
[5] and |6j Between the electrodes we notice the typical structure of a capacitively coupled 
discharge. The plasma is separated in two distinct regions, the plasma bulk and the sheath. 
The maximum of the electron density is approximately 1.9x10^° cm"^. In the sheaths it 
decreases to a much smaller value. A strongly electronegative characteristic can be observed. 
The charged particles are dominated by the oxygen species. The positive ion density is 
almost exclusively due to O^; 0+ and He"*" are neghgible. The negative ions are O3 , O2 and 
O^, in descending order. The electron temperature Te in the plasma bulk is approximately 
2 eV, in the sheath region it is slightly larger due to increased Ohmic heating and possibly 
due to secondaries. 

The effluent of the jet has the character of an afterglow. Both the electron density and 
the electron temperature drop abruptly. This can also be seen from the curves in Fig. [7] 
which show the distribution of the two quantities along the middle axis of the simulation. In 
the jet channel, Ue and are almost constant; after the end of the electrodes (the beginning 
of the effluent) - indicated by the vertical dotted line - both quantities decay exponentially. 
Fig. 8 shows that this feature holds also for all other charged particles except for and O J . 
We conclude from these results that the plasma dynamics in the jet is (with respect to the 
x-coordinate) "local". Generation of plasma is governed by the local electron temperature 
profile which in turn is a functional of the local field distribution; losses are dominated by 
drift and diffusion to the electrodes: The plasma chemistry is so "fast" that transport due 
to the fiow plays a minor role and the solutions become essentially one-dimensional. The 
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agreement of our results with the results performed by Waskoenig et al. by means of a 
ID-fluid model is thus not a coincidence [Hj.) 

The picture changes drastically when we focus on the reactive neutral species. (See Fig. 
9|) For the neutrals the chemistry is "slow", and it takes a considerable time (or distance 
along the flow) for the species to build up. Particularly O, O3 and 02(^A) increase in the jet 
channel. In the effluent, their destruction is also slow; they are transported over a distance 
of a few cm. The density of ozone even increases slightly and remains constant for a certain 
length. The reason for this behavior is that the activated species are not directly generated 
by the plasma electrons but by heavy particle chemistry. A particular case are the ions 
and O3 . The generation of is via electron impact ionization and therefore fast, the 
generation of Oj is slow because it involves ion chemistry. The destruction is slow for both 
species, they are transported about 1 cm into the effluent: and O3 form the end points of 
the positive and negative ion chemistry, respectively, and their direct recombination constant 
is relatively small. 

Further down in the effluent, therefore, only three species play a role: O, O3 and 02(^A). 



Their spatial distribution is shown in Figs. 10, II. and 12 Atomic oxygen reacts with 
molecular oxygen to ozone; its density thus decays exponentially. The rate constant, how- 
ever, is relatively small; 2 cm behind the nozzle of the jet it still has a concentration of 
3x10^^ cm~^. The densities of ozone (2x10^^ cm~^) and excited oxygen (1x10^^ cm~'^) 
remain nearly constant up to the outlet. 

An exhaustive comparison of our simulations with experimental data is currently under- 
way and will be presented in the near future. Nonetheless, some preliminary results may be 



of interest: Fig. 13 shows the profiles of O and O3 along the jet and the effluent, in com- 
parison with experimental results obtained by Ellerweg et al. with quantitative molecular 
beam mass spectrometry and two-photon absorption laser-induced fiuorescence spectroscopy 
[IS] . We find that the simulation results are in excellent agreement with the experiment. 
(Note that we changed the simulation parameters slightly to accommodate the experiment: 
The jet gap size was decreased to 1 mm. The power deposition was kept constant and the 
temperature of the electrodes a very sensitive parameter which controls the production of 
reactive oxygen species in the effluent was set to 330 K.) 
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V. SUMMARY AND CONCLUSION 



In this work we employed the code nonPDPSIM to conduct a numerical simulation 
of an RF-driven atmospheric pressure plasma jet operated with a mixture of helium and 
oxygen. In several aspects our study goes beyond previous work. The two-dimensional 
spatial resolution of the simulation allows to investigate not only the profile of the discharge 
parameters between the electrodes but also their transport through the jet and into the 
environment. We found a pronounced difference in the behavior of the charged particles and 
the neutrals. The charges - e and predominantly O^, 0~, O2 - are governed by a "fast" 
chemistry: In the jet, their production is local in x and their losses are dominated by drift 
and diffusion to the electrodes - this can be understood in a Id-picture which resolves only 
the y-axis [H]. In the effluent, they are virtually absent. The neutrals, in contrast, obey a 
"slow" chemistry. The reactive oxygen species O, O3 and 02(^A) build up along the x-axis 
of the jet and are then flow-transported far into the effluent; their distribution is genuinely 
"two-dimensional". (Actually, it will be "three-dimensional". We believe, however, that 
the restriction of our analysis to two dimensions does not introduce more than quantitative 
errors.) 

In summary, our simulation results indicate that the influence of the jet on the environ- 
ment is primarily of chemical nature. Even in the jet, the density of the activated oxygen 
species is four orders of magnitude higher than the maximum density of the charges, and in 
the effluent the charges vanish completely. This view is in agreement with the experimental 
findings established by various other researchers [T3HT81 [261 l27j . In particular ozone with a 
longer lifetime than atomic oxygen and the very reactive 02(^A) are able to react with sur- 
faces over a long distance [28]. Excited atomic oxygen species, in contrast, plays a relatively 
minor role due to its lower concentration. 

Further work on the /iAPPJ will focus on the transient behavior of the device and also on 
its interaction with a more complex environment. In particular, it is planned to investigate 
the effect of an atmosphere which consists of ambient air (instead of pure helium), and also 
the interaction of the effluent on a material substrate. 
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FIG. 1: Photograph of the micro-atmospheric pressure plasma jet (//APPJ). (Courtesy of Dr. V. 
Schulz-von der Gathen.) The discharge is confined by the top and bottom electrode and the quartz 
glasses at the sides of the jet. The stainless steel electrodes are driven at an RF frequency of 
13.56 MHz and a nominal power of 10 W. The quartz cuvette allows optical access for diagnostic 
purposes. The gas mixture is injected in the jet with a high flow rate (left hand side of the jet), 
while the effluent is formed outside the jet (right hand side of the jet). 
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FIG. 2: Simulation domain with numerical mesh. The triangular unstructured mesh consists of 
approximately 10000 nodes. Around 8000 nodes are located in the plasma region (including the 
effluent), most of them between the top and bottom electrode. The mesh becomes coarser with 
increasing distance to the jet nozzle, above and under the electrodes. On the outer border of 
the simulation domain we employ the Dirichlet boundary condition $ = for the electrostatic 
potential (except at the electrodes). We choose an appropriate distance between the electrodes 
and the boundary to avoid parasitic discharges. 
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FIG. 3: Spatially resolved electron density rie averaged over one RF period in units of cm The 
electron density is confined between the electrodes - it decays rapidly in the effluent. 
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Electron temperature (eV) 
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FIG. 4: Spatially resolved electron temperature Te averaged over one RF period in eV. The electron 
temperature stays almost constant along the jet channel while decreasing exponentially in the 
effluent. 
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FIG. 5: Densities of charged species (total density of positive ions n^^^^i, negative ions n^^^^i and 
electrons e) averaged over one RF period between the electrodes across the jet channel at x = 2 
cm in units of cm~^. The discharge shows a significant electronegative behavior. 
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Electron temperature Iwtween the electrodes (eV) 




FIG. 6: Electron temperature Tg averaged over one RF period between the electrodes across the 
jet channel at x = 2 cm in units of eV. The maxima of the electron temperature are located at the 
sheath edge. 
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Electron density (cm" ) and electron temperature (eV)' 
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FIG. 7: Electron density (in cm~^) and temperature (in eV) in the center between the electrodes 
along the jet channel, averaged over one RF period. The vertical dotted line indicates the edge of 
the electrodes. 
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FIG. 8: Logarithmic-linear plot of densities of charged species along the middle axis of the plasma 
jet. The dotted line at a; = 3 cm indicates the end of the electrodes of the jet. 
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Densities of reactive and metastable species along the jet (cm'^) 
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FIG. 9: Logarithmic-linear plot of the reactive species and metastables densities along the middle 
axis of the plasma jet. The dotted line at x = 3 cm indicates the end of the electrodes of the jet. 
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Atomic oxygen density (cm'^) 
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FIG. 10: Spatial distribution of atomic oxygen density in units of cm~^. After reaching 
maximum of the concentration at the nozzle of the jet the density decays in the effluent. 
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FIG. 11: Spatial distribution of ozone density in units lof cm~^. The ozone density is nearly 
constant across the gap and builds up linearly along the flow. A nearly constant level of the ozone 
concentration maintains in the effluent. 
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FIG. 12: Spatial distribution of metastable molecular oxygen 02(^A) density in units of cm~^. 
The density of metastable oxygen is nearly constant across the gap and grows linearly with respect 
to the lateral coordinate x. In the effluent a slow decay is observed. 
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Densities of reactive oxygen species along the jet (cm' ) 



10^' 



10^' 



10'- 




X position (cm) 



FIG. 13: Comparison of the spatial distribution of the densities of reactive oxygen species in units 
of cm~^ from the simulation with experimental results |18j . 
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